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Global analyses of Parton Distribution Functions (PDFs) have provided incisive constraints on the 
up and down quark components of the proton, but constraining the other flavor degrees of freedom 
is more challenging. Higher-order theory predictions and new data sets have contributed to recent 
improvements. Despite these efforts, the strange quark PDF has a sizable uncertainty, particularly 
in the small x region. We examine the constraints from experiment and theory, and investigate the 
impact of this uncertainty on LHC observables. In particular, we study W/Z production to see how 
the s-quark uncertainty propagates to these observables, and examine the extent to which precise 
measurements at the LHC can provide additional information on the proton flavor structure. 
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I. INTRODUCTION 



A. Motivation 

Partem distribution functions (PDFs) provide the es- 
sential link between the theoretically calculated partonic 
cross-sections, and the experimentally measured physical 
cross-sections involving hadrons and mesons. This link 
is crucial if we are to make incisive tests of the standard 
model, and search for subtle deviations which might sig- 
nal new physics. 

Recent measurements of charm production in neutrino 
deeply- inelastic scattering (DIS), visible as di-muon final 
states, provide important new information on the strange 
quark distribution, s(x), of the nucleon [ll-tlaj. We show 
that despite these recent advances in both the precision 
data and theoretical predictions, the relative uncertainty 
on the heavier flavors remains large. We will focus on the 
strange quark and show the impact of these uncertainties 
on selected LHC processes. 

The production of W/Z bosons is one of the "bench- 
mark" processes used to calibrate our searches for the 
Higgs boson and other "new physics" signals. We will 
examine how the uncertainty of the strange quark PDF 
influences these measurements, and assess how these un- 
certainties might be reduced. 



B. Outline 

The outline of the presentation is as follows. In Sec- 
tion 2, we examine the experimental signatures that con- 
strain the strange quark parton distribution. In Section 3 
we consider the impact of s-quark PDF uncertainties on 
W/Z production at the LHC, and in Section 4 we sum- 
marize our results. Additional details on PDF fits to di- 
muon data at next-to-leading order are provided in the 
Appendix. 



II. CONSTRAINING THE PDF FLAVOR 
COMPONENTS 

A. Extracting the Strange Quark PDF 

In previous global analyses, the predominant informa- 
tion on the strange quark PDF s(x) came from the dif- 
ference of (large) inclusive cross sections for neutral and 
charged current DIS. For example, at leading-order (LO) 
in the parton model one finds that the difference be- 
tween the Neutral Current (NC) and Charged Current 
(CC) DIS i<2 structure function is proportional to the 
strange PDF. Specifically if we neglect the charm PDF 
and isospin-violating terms, we have [17j 
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Figure 1: Relative uncertainty of the strange quark PDF 
as a function of x for Q — 2 GeV. The inner band is 
for the CTEQ6.1 PDF set, and the outer band is for 
the CTEQ6.6 PDF set. The band is computed as the 
envelope of Si(x)/so(x) where Sq(x) is the central PDF 
for each set; for CTEQ6.1, i = [1,40], and for CTEQ6.6, 
i=[l,44]. 



Because the strange distributions are small compared to 
the large up and down PDFs, the s(x) extracted from 
this measurement has large uncertainties. Lacking better 
information, it was commonly assumed the distribution 
was of the form 



s(x) = s(x) ~ k[u(x) + d(x)]/2 



(2) 



with 
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[s(x) + s(x)} 
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This approach was used, for example, in the CTEQ6.1 
PDFs [13| • In Figure Q] we show the relative uncertainty 
band of the strange quark PDF for the 40 CTEQ6.1 PDF 
error sets relative to the central value. We observe that 
over much of the x-range the relative uncertainty on the 
strange PDF is < 5%. The relation of Eq. © tells us that 
this uncertainty band in fact reflects the uncertainty on 
the up and down sea which is well constrained by DIS 
measurements; this does not reflect the true uncertainty 
of s(x). 

Beginning with CTEQ6.6 PDFs (Hi the neutrino- 
nucleon dimuon data was included in the global fits to 
more directly constrain the strange quark; thus, Eq. @ 
was not used, and two additional fitting parameters were 
introduced to allow the strange quark to vary indepen- 
dently of the up and down sea. We also display the rela- 
tive uncertainty band for the CTEQ6.6 PDF set in Fig.Q] 
We now observe that the relative error on the strange 
quark is much larger than for the CTEQ6.1 set, partic- 
ularly for x < 0.01 where the neutrino-nucleon dimuon 
data do not provide any constraints. We expect this is a 
more accurate representation of the true uncertainty. 

This general behavior is also exhibited in other global 
PDF sets with errors (MSI. For example, the NNPDF 
collaboration uses a parameterization-free method for ex- 
tracting the PDFs; they observe a large increase in the 



s(x) uncertainty in the small x region which is beyond the 
constraints of the i^-DIS experiments. (Cf. in particular 
Fig. 13 of Ref (p.) 

Thus, there is general agreement that the strange 
quark PDF is poorly constrained, particularly in the 
small x region. 



B. Constraints from CCFR and NuTeV 

The primary source of information on the strange 
quark at present comes from high-statistics neutrino- 
nucleon DIS measurements; in particular, the CCFR and 
NuTeV dimuon experiments have been used to deter- 
mine the stran ge quark PDF with improved accuracy 
Is IE S Ha, l24l - [26fl . Neutrino induced dimuon produc- 
tion (yN — > f.i + [J,~X) proceeds primarily through the 
Cabibbo favored s — > c or s — > c subprocess. Hence, 
this provides information on s and s directly; this is in 
contrast to AF 2 of Eq. (fT). CCFR has 5030 v and 1060 
v di-muon events, and NuTeV has 5012 v and 1458 v 
di-muon events, and these cover the approximate range 
x ~ [0.01,0.4]. Additionally, NuTeV used a sign-selected 
beam to separate the v and v events in order to sepa- 
rately extract s(x) and s(x). 



1. Constraints on s + s 

In Table Q] we illustrate how the bulk of the data 
used in the global fits are relatively insensitive to the 
strange quark distribution. The first column (labeled 
"CTEQ6M") lists the X 2 /DoF for a variety of data sets 
used in the CTEQ6M fit Q2. We have also shown the 
CCFR and NuTeV dimuon data sets in the table, but 
these were not used in the CTEQ6M fit. The second col- 
umn (labeled "Free") lists results of refitting all the data 
- including the dimuon data - with a flexible strange 
quark PDF instead of imposing the relation of Eq. ([2); 
this allows the strange quark PDF to accommodate the 
dimuon data. Comparing the two columns, we observe 
that the change of the strange PDF allowed for a greatly 
improved fit of the dimuon data, while the other data 
sets are virtually insensitive to this change |j 

This exercise demonstrates that most of the data sets 
of the global analysis are insensitive to the details of the 
strange quark PDF. 



X 2 /DoF 


CTEQ6M 


Free 


CCFR v dimuon 


1.02 


0.72 


CCFR v dimuon 


0.58 


0.59 


NuTeV v dimuon 


1.81 


1.44 


NuTeV v dimuon 


1.48 


1.13 


BCDMS F,f 


1.11 


1.11 


BCDMS F$ 


1.10 


1.11 


HI 96/97 


0.94 


0.94 


HI 98/99 


1.02 


1.03 


ZEUS 96/97 


1.14 


1.15 


NMC Fl 


1.52 


1.49 


NMC Fi/Fl 


0.91 


0.91 


CCFR F 2 


1.70 


1.88 


CCFR F3 


0.42 


0.42 


E605 


0.82 


0.83 


NA51 


0.62 


0.52 


CDF £ Asym 


0.82 


0.82 


E866 


0.39 


0.38 


DO Jets 


0.71 


0.67 


CDF Jets 


1.48 


1.47 


TOTAL x 2 


2173 


2133 



Table I: We display the x 2 / DoF for selected data sets 

using the CTEQ6M PDF set [Hj|, and a variant of this 

(labeled "Free") which allows for a modified strange quark 

PDF to accommodate the neutrino dimuon data. 



the relative x 2 for the dimuon and "Inclusive- 1" data sets 
as a function of the strange asymmetry [S~\, where 



\S~ 



x[s(x) — s(x)] dx 



(3) 



The "Inclusive-I" data sets (cf., Ref. [27|]) contain the data 
that is sensitive to [S~]; specifically, the data sets are 
a) the neutrino xFs data from CCFR and CDHSW as 
this is proportional to the difference of quark and anti- 
quark PDFs, and b) the CDF VF-asymmetry measure- 
ment which can receive contributions from the sg — > Wc 
subprocess. Figure [2] clearly shows that the dimuon data 
provides the strongest constraints on the strange asym- 
metry [S~]. 



2. Constraints on s 



HERMES 



The dimuon data can also provide information on the 
s(x) and s(x) quark PDFs separately. In Fig.[5]we display 



1 The one exception is the CCFR F2 which is mildly sensitive to 
the strange quark PDF via Eq. JTJ. 



The HERMES experiment measured the strange 
PDF via charged kaon production in positron-deuteron 
DIS [23|, these results are displayed in Fig. [3j For com- 
parison, the strange quark and total sea distributions 
from CTEQ6L are also plotted. 

The HERMES data suggests that the x-dependence of 
the strange quark distribution is quite different from the 




[S"]xl0 4 

Figure 2: We plot X /Xo f° r the dimuon and the 
"Inclusive-I" data sets evaluated as a function of the 
strange asymmetry [S~] x 10 4 . The fits are denoted with 
■ for the dimuons and A for "Inclusive-I". Quadratic ap- 
proximations to the fits are displayed by the solid (red) 
line for the dimuons and the dashed (green) line for 
"Inclusive-I". 



form assumed for the CTEQ6 set. In particular, they 
obtain a strange quark distribution that is suppressed in 
the region x > 0.1 but then grows quickly for x < 0.1 
and exceeds the CTEQ6L value in the small x region by 
more than a factor of two. 

To gauge the compatibility of this result with the dis- 
played PDFs, we can replace the initial s{x) distribution 
with the form preferred by HERMES, and then evalu- 
ate the shift of the x 2 with this additional constraint. A 
preliminary investigation with this procedure indicates 
that the HERMES s(x) distribution could strongly influ- 
ence two data sets of the global fits. The first set is the 
neutrino- nucleon dimuon data which controls s(x) in the 
intermediate x region. The second set is the HERA mea- 
surement of i<2 in the small x region where the statistical 
errors are particularly small. 

In Fig. [3] we also show xS(x) from CTEQ6.6; while 
the HERMES data are below the CTEQ6.6 result in the 
x ~ 0.1 region, they agree quite well at both the higher 
and lower x values. 

While these comparisons are sufficient to gauge the 
general influence of the Hermes result, a complete anal- 
ysis that includes the Hermes data dynamically in the 
global fit is required to draw quantitative conclusions. 
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Figure 3: The strange parton distribution xS(x) = 
x[s(x) + s(x)] from the measured Hermes multiplicity 
for charged kaons evolved to Q 2 = 2.5 GeV 2 . The 
solid green curve is a Hermes 3-parameter fit: S(x) — 
x -o.924 e -a:/o.o404^_ a ,^ y^ hashed \)\ ue curve is the sum 

of light anti-quarks x(u+d) from CTEQ6L, the blue solid 

curve is xS(x) from CTEQ6L, and the red solid curve is 

the xS(x) from CTEQ6.6. Hermes data points and fit 

are from Ref. (H. 



D. CHORUS 



The CHORUS experiment [29(43 1| measured the neu- 
trino structure functions F2, xF^, R in collisions of sign 
selected neutrinos and anti-neutrinos with a lead target 
(lead-scintillator CHORUS calorimeter) in the CERN 
SPS neutrino beamline. They collected over 3M v^ and 
1M £/p charged current events in the kinematic range 
0.01 < x < 0.7, 0.05 < y < 0.95, 10 < E v < 100. 

This data was analyzed in the context of a global fit 
in Ref. [H which was based on the CTEQ6.1 PDFs. 
This analysis made use of the correlated systematic errors 
and found that the CHORUS data is generally compati- 
ble with the other data sets, including the NuTeV data. 
Thus, the CHORUS data is consistent with the strange 
distribution extracted in CTEQ6.1. 



E. NOMAD 

The NOMAD experiment measured neutrino-induced 
charm dimuon production to directly probe the s-quark 
PDF [33l435|]. Protons from the CERN SPS synchrotron 
(450 GeV) struck a beryllium target to produce a neu- 
trino beam with a mean energy of 27 GeV. NOMAD used 
an iron-scintillator hadronic calorimeter to collect a very 
high statistics (15K) neutrino-induced charm dimuon 
sample [34} . 

Using kinematic cuts of E^i^E^ > 4.5 GeV, 15 < 
E v < 300 GeV, and Q 2 > 1 GeV 2 NOMAD performed 
a leading-order QCD analysis of 2714 neutrino- and 115 



anti-neutrino-induced opposite sign dimuon events [33| . 
The ratio of the strange to non-strange sea in the nucleon 
was measured to be k = 0.48^ ' 07 ^ 012 ; this is consistent 
with the values used in the global fits, c.f., Fig. [5] 

The data analysis is continuing, and it will be very 
interesting to include this data set into the global fits as 
the large dimuon statistics have the potential to strongly 
influence the extracted PDFs. 



F. MINERz/A 

The cross sections in neutrino DIS experiments from 
NuTeV, CCFR, CHORUS and NOMAD have been mea- 
sured using heavy nuclear targets. In order to use these 
measurements in a global analysis of proton PDFs, these 
data must be converted to the corresponding proton or 
isoscalar results |36T - t42| . For example, the nuclear cor- 
rection factors used in the CTEQ6 global analysis were 
extracted from £^N DIS processes on a variety of nuclei, 
and then applied to uN DIS on heavy nuclear targets. 
In a series of recent studies it was found that the £^N 
nuclear correction factors could differ substantially from 
the optimal vN nuclear correction factors [39h43| . 

Furthermore, the nuclear corrections depend to a cer- 
tain degree on the specific observable as they contain 
different combinations of the partons; the nuclear correc- 
tion factors for dimuon production will not be exactly 
the same as the ones for the structure function F-x or 
i*3 . The impact of varying the nuclear corrections on the 
strange quark PDF has to be done in the context of a 
global analysis which we leave for a future study. 

The MINERVA experiment has the opportunity to 
help resolve some of these important questions as it can 
measure the neutrino DIS cross sections on a variety of 
light and heavy targets. It uses the NuMI beamline at 
Fermilab to measure low energy neutrino interactions to 
study neutrino oscillations and also the strong dynamics 
of the neutrino-nucleon interactions. MINERVA com- 
pleted construction in 2010, and they have begun data 
collection. MINERVA can measure neutrino interactions 
on a variety of targets including plastic, helium, carbon, 
water, iron, and lead. For 4 * 10 20 Protons on Target 
(POT) they can generate over 1M charged current events 
on plastic. 

These high statistics data on a variety of nuclear tar- 
gets could allow us to accurately characterize the nuclear 
correction factors as a function of the nuclear A from 
helium to lead. This data will be very useful in resolv- 
ing questions about the nuclear corrections, and we look 
forward to the results in the near future. 



G. CDF & DO 

At the Tevatron, the CDF 0] and DO [45] collabora- 
tions measured Wc final states in pp at ^/S = 1.96 TeV 
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Figure 4: k(x,Q) vs. x for Q = 1.5 GeV for a se- 
lection of PDFs, where k(x,Q) is defined in Eq. (QJ. 
The curves (top to bottom) are CTEQ6.6 (solid, red), 
CTEQ6.5 (dotted, black) and CTEQ6.1 (dashed, pur- 
ple). The wider (blue) band represents the uncertainty 
for CTEQ6.6 as computed by Eq. ([S]), the inner (green) 
band represents uncertainty given by the envelope of 
k(x, Q) values obtained with the 44 CTEQ6.6 error sets. 



using the semileptonic decay of the charm and the corre- 
lation between the charge of the W and the charm decay. 
Additionally, a recent study has investigated the im pac t 
of the IF+dijet cross section on the strange PDF |4a |. 
These measurement are especially valuable for two rea- 
sons. First, there are no nuclear correction factors as the 
initial state is p or p. Second, this is in a very different 
kinematic region as compared to the fixed-target neutrino 
experiments. Thus, these have the potential to constrain 
the strange quark PDF in a manner complementary to 
the uN DIS measurements; however, the hadron-hadron 
initial state is challenging. Using approximately 1 /6 _1 
of data, both CDF and DO find their measurements to be 
in agreement with theoretical expectations of the Stan- 
dard Model. Updated analyses with larger data sets are 
in progress and it will be interesting to see the impact of 
these improved constraints on the strange quark PDF. 



H. Strange Quark Uncertainty 

The combination of the above results underscores the 
observation that our knowledge of the strange quark is 
limited. To illustrate this point in another manner, in 
Fig. [5] we display k(x, Q) for a selection of PDF sets. 
Here, we define 



k(x,Q) 



s{x 7 Q) 



[u(x,Q) + d(x,Q)]/2 



(4) 



which is essentially a differential version of the k param- 
eter of Eq. ([2]); this allows us to gauge the amount of the 
strange PDF inside the proton compared to the average 
up and down sea-quark PDFs. If we had exact SU(3) 




v^ 



0.0001 



Figure 5: k(x, Q) vs. x showing the evolution from low to 
high scales. The solid (red) lines are for CTEQ6.6, and 
the dashed (purple) lines are for CTEQ6.1. The lower 
pair of lines (red and purple) are for Q = 1.5 GeV and 
the upper for Q = 80 GeV. 



symmetry we would expect u = d = s and k(x, Q) ~ 1. 
As the strange quark is heavier than the up and down 
quarks, we expect this component to be suppressed rel- 
ative to the up and down quarks, and we would predict 
u ~ d > s which would yield k(x, Q) < 1. Thus, k(x, Q) 
is a measure of the SU(3) breaking across the x and Q 
range. 

In Fig. H we observe that the CTEQ6.1 and CTEQ6.5 
PDF sets have k(x, Q) ~ 1 /2\ this was by design as the 
constraint of Eq. (|5J) was used to set the initial s(x) distri- 
bution. The exception is CTEQ6.6 which did not impose 
Eq. I2J); we observe that this set has k(x,Q) ~ X J2 for 
x ~ 0.1 (where the dimuon DIS data has smaller un- 
certainties) , but is a factor of two larger than the other 
PDF sets for small x values. In Fig. [H we also show the 
uncertainty on s(x) computed as [12| 



AX = - 
2 



\ 



A> 



*£[x(st)-x(s-)}- 



(- r >) 



which is shown as a (blue) bandjj this results in a band 
which is larger than simply taking the spread of the 44 
CTEQ6.6 error PDFs (green band). 

In order to show the effect of the DGLAP evolu- 
tion on the strange distribution, we display k(x, Q) for 
CTEQ6.1 and CTEQ6.6 at both a low and high Q scale 
in Fig. [5J As we explore higher scales, the production of 
s(x) by gluon splitting moves k(x, Q) toward the SU(3)- 
symmetric limit. This trend is especially pronounced at 
low x-values. Thus, as the LHC W/Z production is cen- 
tered in the range x ~ 0.01, we will be particularly inter- 
ested in the k(x, Q) changes in this region. 



2 In Eq. J5J, X is the observable, S. are the error PDF sets for 
eigenvalue i, and N p is the number of eigenvalues. For CTEQ6.5 



S =14TeV 7TeV 
1 



0.1 



1.96TeV 



0.01 
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N„ 



20, and for CTEQ6.6 N p = 22. 



10~ 4 0.001 0.01 0.1 
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Figure 6: Parton momentum fractions xa and xb ac- 
cessible in W and Z boson production in the Tevatron 
Run-2 (VS = 1.96 TeV), and at the LHC (VS = {7, 14} 
TeV) . The accessible ranges of xa and Xb are shown by 
the solid lines. The contours of the constant rapidity y 
are shown by the inclined dotted lines. 



These results reflect the relevant arrange of the con- 
straints on the strange quark PDF, and how they depend 
on the Q-scale. In the next section we will investigate the 
implications of this uncertainty on the Drell-Yan W/Z 
boson production at the LHC. 



III. IMPLICATIONS FOR DRELL-YAN W/Z 
PRODUCTION AT THE LHC 



The Drell-Yan production of W ± and Z bosons at 
hadron colliders can provide precise measurements for 
electroweak observables such as the W boson mass |47l . 
|48| and width, the weak mixing angle in 7* jZ produc- 
tion [49j, and the lepton asymmetry in W production. 
These results can measure fundamental parameters of 
the Standard Model (SM) and constrain the Higgs boson 
mass. If a Higgs boson is found at the LHC, Drell-Yan 
W/Z boson production will help in the search for devia- 
tions of the SM and to reveal new physics signals |50l— f52j| . 
For instance, new heavy gauge bosons could be discov- 
ered in the invariant lepton distribution or new particles 
and interactions might leave a footprint in the Peskin- 
Takeuchi S and T parameters |53j |. 

Furthermore, the W and Z boson cross section "bench- 
mark" processes are intended to be used for detector 
calibration and luminosity monitoring [54|; to perform 
these tasks it is essential that we know the impact of the 
PDF uncertainties on these measurements. The impact 
on these benchmark processes, and the Higgs boson pro- 
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(a) da/dy for W (left), W + (middle), Z° (right) boson production at the Tevatron. 
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(b) rfo-/di/ for W~ (left), VK+ (middle), Z° (right) boson production at the LHC with Vs = 7 TeV. 
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(c) dcr/dy for W - (left), Vy+ (middle), Z° (right) boson production at the LHC with </S = 14 TeV. 

Figure 7: Partonic contributions to the differential cross section of on-shell W/Z boson production at LO as a 

function of the vector boson rapidity. Partonic contributions containing a strange or anti-strange quark are denoted 

by (red) dashed and (blue) dot-dashed lines. The solid lines show the total contribution. 



duction, were studied in Refs. |55| - |57| . In the following, 
we will investigate the influence of the PDFs on the ra- 
pidity distributions of the Drell-Yan production process. 
Conversely, it may be possible to use the W/Z production 
process to further constrain the parton distribution func- 
tions in general, and the strange quark PDF in particular. 
As noted in Ref. |49|, when looking for new physics sig- 
nals it is important not to mix the information used to 
constrain the PDFs and the new physics as this would 
lead to circular reasoning. 

As we move from the Tevatron to the LHC scatter- 
ing processes, the kinematics of the incoming partons 
changes considerably; in Fig. [6] we show the momentum 
fractions xa and xb of the incoming parton A and par- 
ton B for the Tevatron Run-2 (\/S = 1.96 TeV) and 
the LHC with VS = 7 TeV and s/S = 14 TeV. The 
solid (red) lines show the range of xa and Xb probed by 
W^- and Z boson production. At the Tevatron, values 
of %a,b down to 2 x 10 -3 are probed for large rapidities 
of yw/z = 3. However, at the LHC much smaller values 



of xa and xb become important due to the larger CMS 
energy and broader rapidity span. For vS = 7 TeV, the 
PDFs are probed for x-values as small as 2 x IO -4 for 
rapidities up to ~ 4.5. With \/S = 14 TeV, even larger 
rapidities of y ~ 5 and smaller values of Xa/b of 4 x 10~ 5 
might be reached. 



A. LHC Measurements 

The importance of the PDF uncertainties to the LHC 
measurements was already evident in the 2010 and pre- 
liminary 2011 data. 

ATLAS presented measurements of the Drell-Yan W/Z 
production at the \/S = 7 TeV with 35 pb- 1 [58]. These 
results include not only the measurement of total cross 
section and transverse distributions, but also a first mea- 
surement of the rapidity distributions for Z — > l + l~ as 
well as W + —> l + vi and W~ — > l~v\. Additionally, 
ATLAS has used W/Z production to infer constraints 



s 



on the strange quark distribution, and they measure 



-,+0.25 
-0.28 



at Q z 



1.9 GeV^ and 



r s = 0.5(s + s)/d = LOO] 
x = 0.023 (Ei|. 

CMS has measured the rapidity and transverse mo- 
mentum distributions for Z — > l + l~ production [60] and 
inclusive W/Z production [6l| using 36 pb~ x of data. 
Additionally CMS has measured the weak mixing an- 
gle [49J, the forward-backward asymmetry in j*/Z pro- 
duction I62| , and the lepton charge asymmetry in W pro- 
duction [6l [H| . 

LHCb has measured the W charge asymmetry in 
Refs. J65| and [66| • These measurements show, already 
with these data samples, the PDF uncertainties are im- 
portant and can be the leading source of measurement 
uncertainty. 

Additionally, CMS has analyzed W + c production 
which is directly sensitive to the s and s contribution 
of the proton; the results for the 36 p6 -1 data sample are 
given in Ref. |67J . 



B. Strange Contribution to W/Z Production 

Because the proton-proton LHC has a different initial 
state and a higher CMS energy than the Tevatron, the 
relative contributions of the partonic subprocesses of the 
W^ jZ production change significantly. At the LHC, the 
contributions of the second generation quarks {s,c} are 
greatly enhanced. Additionally, the W + and W~ rapid- 
ity distributions are no longer related by a simple y — > — y 
reflection symmetry due to the pp initial state. In Fig- 
ure [7] we display the contributions from the different par- 
tonic cross sections which contribute to W^ and Z pro- 
duction at LO. 

Figure [7a| shows the rapidity distribution at the Teva- 
tron. For W + (W~) production, the ud (du) channel 
(dotted black lines) contributes 90% of the cross section, 
while in Z production the uu (dotted black line) and 
dd (dot-dashed black line) subprocesses contribute 93% 
of the cross section. The first generation quarks {u, d} 
therefore dominate the production process while contri- 
butions from strange quarks (red dashed and blue dot- 
dashed lines) are comparably small with 9% for W and 
5% (ss) for Z boson production. 

At the LHC, subprocesses containing strange quarks 
are considerably more important as shown in Fig. I7bl for 
a CMS energy of 7 TeV and in Fig. [7c] for 14 TeV. For 
W~ production (left plots), the (blue) dot-dashed lines 
show the su channel while the (red) dashed lines show 
the sc contribution. At 14 TeV the sc —J- W~ subprocess 
contributes 28% to the cross section, while the su — > W~ 
subprocess contributes only 2% as this is suppressed by 
the off-diagonal CKM matrix entry. For W + produc- 
tion channels, the su channel (blue dot-dashed lines) con- 
tributes only 2%, while the sc channel (red dashed lines) 
yields 21%. Notice, the absolute value of the sc — > W~ 
and sc — > W + contributions are the same, however the 
relative contribution is smaller for W + production due to 



the larger up-quark valence contribution in the ud —> W + 
subprocess as compared to ud — > W~ . 

The rapidity distributions of the total W~ and W + 
boson production differ markedly at the LHC because of 
the different valence quark contributions from u and d. 
This effect is also present in the su —$■ W~ and su — > W + 
(blue dot-dashed lines) subprocess. We will comment 
more on this feature in the following subsection. 

Comparing Figs. [7a| with Figs. [7cl we note the LHC 
explores a much larger rapidity range. For channels con- 
taining strange quarks, \yw/z\ can be measured up to 
y ss 4.5 at the LHC, compared to y « 2.5 at the Teva- 
tron; therefore smaller values of x of the strange quark 
distribution can be probed. 

Additionally, as at the Tevatron, we can use W~ pro- 
duction to probe the strange quark PDF while using W + 
production to probe the anti-strange PDF. 

While the LO illustration of Fig. [7J provides a useful 
guide, in Fig. [5] we display the strange quark contribution 
to the differential cross section d 2 a/dM/dy of on-shell 
W~ , W + , Z boson production computed at NNLO using 
the VRAP program [68( . We display the LHC results for 
W ± and Z with y/S of both 7 TeV and 14 TeV, where 
the (yellow) band represents the strange-quark initiated 
contributions to the total differential cross section. 

The figures impressively highlight the large contribu- 
tion of the strange and anti-strange quark subprocesses 
at the LHC. Consequently it is essential to constrain the 
strange PDF if we are to make accurate predictions and 
to perform precision measurements. Figure[5]also demon- 
strates clearly the very different rapidity profiles of the 
strange quark (arising from the sea-distribution) com- 
pared to the u and d quark terms which are dominated 
by the valence distributions. This property is most ev- 
ident for the case of W + production. Here, the domi- 
nant ud contribution has a twin-peak structure due to 
the harder valence distribution, while the cs distribution 
has a single-peak centered at y = 0. The total distribu- 
tion is then a linear combination of the twin-peak and 
single-peak distributions, and these are weighted by the 
corresponding PDF. 

Therefore, a detailed measurement of the rapidity dis- 
tribution of the W^/Z bosons can yield information 
about the contributions of the s quark relative to the u, d 
quarks. As this is a relative measurement, rather than 
an absolute cross section measurement, it is reasonable 
to expect that this could be achieved with high precision 
once sufficient statistics are collected. Consequently, this 
is an ideal measurement where the LHC data could lead 
to stronger constraints on the PDFs. 



PDF Uncertainty of the W/Z rapidity 
distributions 



To estimate the influence of the PDF uncertainties 
(and in particular the strange quark PDF) on the W/Z 
production process and its differential distributions at the 
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(a) d?a/dM/dy in pb/GeV for pp->W~ +X (left), pp -> W+ + X (middle) and pp -> Z, 7* + X (right) production at the LHC 
for 7 TeV with CTEQ6.6 using the VRAP program J6J| at NNLO (for M = m w and M = m z respectively). 
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(b) d 2 a/dM/dy in pb/GeV for pp -> W~ + A" (left), pp -> VK+ +X (middle) and pp -> Z,7* + X (right) production at the LHC 
for 14 TeV with CTEQ6.6 using the VRAP program [gjj at NNLO (for M = m w and M = m z respectively). 

Figure 8: Contribution of the strange quark to W^ JZ production at the LHC. 



LHC, we will use the different PDF sets within CTEQ6.6 
as well as compare the sets of different PDF groups. 

In Figure [Ha] we display the differential cross section 
d 2 a/dM/dy for W ± /Z boson production at the LHC at 
\/~S = 7 TeV using the 44 error PDF sets of CTEQ6.6. 
To better resolve these PDF uncertainties, we plot the 
ratio of the differential cross section d 2 a / dM / dy com- 
pared to the central value in Fig. I9bl We observe that 
the uncertainty due to the PDFs as measured by this 
band is between ±3% and ±4% for central boson rapidi- 
ties of — 3 < yw/z < +3- For larger rapidities, the PDF 
uncertainties increase dramatically, but the cross section 
vanishes. 

For comparison, in Fig. [5c] we display the (yellow) band 
of CTEQ6.6 error PDFs together with the results using 
other contemporary PDF sets. The (yellow) band shows 
the span of the 44 CTEQ6.6 error PDFs of Fig. [90 and 
the solid lines show the rapidity distribution from the 
selection of PDFs; all have been scaled to the central 
value for the CETQ6.6 set[j We observe that the choice 
of PDF sets can result in differences ranging up to ± 8 % 
for —2 < yw/z < +2 and even up to ±10% for —3 < 
\yw/z\ < 3, which is well beyond the ±3% and ±4% 



3 Here, we are more interested in the general span of these dif- 
ferent PDFs rather than the specific sets and values. For 
reference the specific curves are: MSTW2008 l2(j| (magenta) . 
NNPDF El (blue), ABKM09 01 (gray), CT10[nl (purpleV 
CTEQ6.5 I7U (black), CTEQ6.1 lial (green), HERAPDF10 Ym 
(orange), MRST2004 [ll (red). 



range displayed in Fig. I9b| note the different scales used 
in Fig. [9b] and Fig. [9c] However, if we compute the PDF 
uncertainty band using Eq. $5§ as specified by Ref. (l2| 
we find an estimated uncertainty of ~ 15% (depending on 
the rapidity) which generally does encompass the range 
of PDFs displayed in Fig. [5c] 

While the band of error PDFs provides an efficient 
method to quantify the uncertainty, the range spanned 
by the different PDF sets illustrates there are other im- 
portant factors which must be considered to encompass 
the full range of possibilities. 



D. Correlations of the W/Z rapidity distributions 

The leptonic decay modes of the W/Z bosons provide a 
powerful tool for precision measurements of electroweak 
parameters such as the W boson mass. As the leptonic 
decay of the W boson contains a neutrino (W — > tv), 
this process must be modeled to account for the missing 
neutrino. The W mass can then be measured by study- 
ing the transverse momentum distribution of the decay 
lepton £ or the transverse mass of the £v pair. Perform- 
ing this measurement, the Drell-Yan Z boson production 
process is used to calibrate the leptonic W process be- 
cause the Z can decay into two visible leptons Z — > £ + £~ . 
This method works to the extent that the production pro- 
cesses of the W and Z bosons are correlated. 

One possible measure to gauge the correlation of the 
PDF uncertainty is the ratio of the W and Z boson dif- 



10 




(a) d 2 a/dM/dy in pb/GeV for pp -> VK~ + X (left), pp -> W+ +X (middle) and pp -> Z,7* +X (right) production at the LHC 

for 7 TeV with CTEQ6.6 using the VRAP program jgjj at NNLO. 
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(b) d 2 a/dM/dy for on-shell {W - , W+, Z} production at the LHC for v^ = 7 TeV with CTEQ6.6 using the VRAP program at 

NNLO, scaled by the central value. 





2 4 -4-2 

y y 

(c) d 2 a/dM/dy for {W~ , W+, Z} production at the LHC for %/S = 7 TeV with a selection of PDFs using the VRAP program at 

NNLO. The (yellow) band is for the CTEQ6.6 set |1S| . and the other curves are for the central values of different PDF sets (see 

text). All plots are scaled by the central value for the CETQ6.6 set. Note the scale of this figure is larger than for Fig. I9bl 

Figure 9: PDF uncertainty bands for on-shell W~(left plots) , W + (middle plots) and Z(right plots) production at 

the LHC for VS = 7 TeV. 



ferential cross section. We compute d 2 a / dM / dy for W ± 
compared to Z, and divide by the central PDF results to 
see the uncertainty band on a relative scale. Schemati- 
cally we define: 



R ± = 



daiW^ 
da(Z) 



I 



daiW^ 



da{Z) 



(6) 



where the "0" subscript denotes the "central" PDF set. 
The resulting distributions are displayed in Fig. HOal for 
W~ production and in Fig. llObl for W + production. The 
left plot in each figure shows the distributions for the 
CTEQ6.5 PDF set, and the right plots the distributions 
for CTEQ6.6. We observe that the uncertainty band is 
generally ±1% for central rapidities of —2 < yw/z < +2; 
this is smaller than in the previous case, where the ab- 
solute uncertainty was investigated. For larger rapid- 



ity ( \yw/z\ > 2) the uncertainty band exceeds the ±1% 
range of the plot. 

In Fig. llOcl we plot the sum of the differential W + 
and W~ cross sections with respect to the differential Z 
boson cross section, again normalized to the distribution 
of the central PDF set. We define: 



R 



da{W+ + W 
da{Z) 



/ 



da{W+ + W- 
da{Z) 



(7) 



for both the CTEQ6.5 and CTEQ6.6 PDFs. 

The contrast in Fig.HUclis striking. For the CTEQ6.5 
PDFs, we observe that VF ± and Z processes are strongly 
correlated, while for the CTEQ6.6 the spread of the PDF 
band is substantially larger. For example, the double 
ratio for CTEQ6.5 has a spread of approximately ±0.2% 
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(a) Double ratio Rr as defined in Eq. JBJ for the LHC for y/S = 7 TeV with CTEQ6.5 
(left) and CTEQ6.6 (right) at NNLO. 
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(b) Double ratio R+ as defined in Eq. (6) for the LHC for y/S = 7 TeV with CTEQ6.5 
(left) and CTEQ6.6 (right) at NNLO. 
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(c) Double ratio R as defined in Eq. Q for the LHC for y/S = 7 TeV with CTEQ6.5 
(left) and CTEQ6.6 (right) at NNLO. 

Figure 10: Ratios of the differential W ± and Z production cross section as denned in Eq. ([6]) and Eq. at the LHC 

forVS = 7 TeV. 



within the central rapidity range of —3 < yw/z < +3, 
while the uncertainty for CTEQ6.6 is much wider in this 
rapidity region. 

The primary difference that is driving this result is the 
different strange PDF. For CTEQ6.5 the strange quark 
was defined by Eq. ([2]) while CTEQ6.6 introduced two 
extra fitting parameters which allowed the strange PDF 
to vary independently from the up and down sea. Thus, 
the uncertainty of the CTEQ6.6 distributions more ac- 
curately reflects the true uncertainty. 

Another means to see how the additional freedom of 
the strange quark introduces a decorrelation of the W 
and Z processes is evident in Fig. [TT] which displays the 
correlation of the W ± and Z boson cross sections for a 
selection of CTEQ PDFs. Except for CTEQ6.6, all the 
PDFs make use of Eq. © and yield results that lie along 
a straight line in the {cr\y, &z} plane. Because CTEQ6.6 



does not use Eq. ©, the freedom of the strange quark 
PDF is reflected in the freedom of the W ± and Z cross 
sections values. 

The above examples demonstrate the subtle features 
inherent in evaluating the PDF uncertainties. For pre- 
cision measurements it is important to better constrain 
the parton distributions at the LHC, in particular the 
strange and anti-strange quark PDFs. 



IV. CONCLUSION 

We have investigated the constraints of the strange 
and anti-strange PDFs and their impact on the Drell-Yan 
W/Z boson production at the LHC. 

Specifically, we observe that the strange quark is rather 
poorly constrained, particularly in the low x region which 
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Figure 11: Correlation of the W^ and Z cross sections 
for a selection of PDF sets. Figure taken from Ref. fl9l | . 



is sensitive to W/Z production at the LHC. Improved 
analyses from neutrino DIS measurements could help re- 
duce this uncertainty. Conversely, precision measure- 
ments of W/Z production at the LHC may provide input 
to the global PDF analyses which could further constrain 
these distributions. 

In particular, the rapidity distribution of the W/Z 
bosons provides an incisive measure of the mix of va- 
lence and sea quarks, and the prospect of measuring this 
at the LHC in the near future is excellent. 



APPENDIX 

A. PDF Fits with the Dimuon Data 

We have repeated the LO analysis of Ref. [27j and ex- 
tended this using the NLO calculation for dimuon pro- 
duction [74l . l7a | . We have performed a series of fits to the 
data which includes the dimuon data. The results of the 
{A, B + , B, B~ , C} fitQ using the LO dimuon analysis are 
shown in Table [TTl and those with the NLO analysis are 
in Table IIII1 The fits are sorted left-to-right by the in- 
tegrated strange-quark asymmetry [S~] (scaled by 10 4 ). 
The cells display the \ 2 relative to Xo for the indicated 
data set where we choose Xo to be the % 2 value from 
the LO-.B fit; this allows us to compare the incremental 
changes as we shift [S~] and alter the constraints. The 
values in parentheses are the x 2 /DoF for each data sub- 
set. The B-fit is the overall best fit to the data. The 



4 We follow the methodology and notation of Ref. [27]] • See this 
reference for details. 



B + and B~ sets modify the B fit using the Lagrange 
multiplier method to determine the ranges of the [S~] 
parameter defined in Eq. ((3|). 

For example, the LO B + fit demonstrates that we can 
increase [S~] from 15.98 to 54.85, but the dimuon \ 2 
increases by 33/174 ~ 21% while the overall \ 2 increases 
by only 39/2465 ~ 2%; thus, the shift of [S~] is strongly 
constrained by the dimuon data, and the remaining data 
are relatively insensitive to this quantity. 

Comparing the NLO-B fit to the LO-i? fit we note that 
X 2 /DoF has decreased both for the dimuon set and the 
entire data set; while this decrease is not dramatic, it 
is encouraging to see that the proper NLO treatment of 
the data results in an improved fit. As before, we observe 
that for the NLO B + fit, we can increase [S~] from 13.72 
to 63.75, but the dimuon x 2 increases by 92/174 ~ 53% 
while the overall % 2 increases by only 106/2465 ~ 4%; 
again, the shift of [S~] is primarily constrained by the 
dimuon data. 

In Figure [2] we have plotted the ratio x 2 /Xo_f° r the 
individual dimuon and "Inclusive-I" data sets |27| evalu- 
ated for a series of NLO S-fits as a function of the strange 
asymmetry [S~] x 10 4 . This plot allows us to see the con- 
tribution of each data set as we shift the strange asym- 
metry. Again the "Inclusive-I" data sets are essentially 
unchanged as the treatment of the dimuons only affects 
these data indirectly. As before, this data set is mildly 
sensitive to the dimuons, and weakly prefers larger values 
of [S-]. 
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LO Dimuon 


#pts 


B+ 


A 


B 


C 


B~ 


[S~] x 10 4 


- 


54.85 


31.18 


15.98 


10.32 


-17.72 


Entire Data 


2465 


1.015 (1.100) 


1.001 (1.085) 


1.000 (1.084) 


1.002 (1.086) 


1.018 (1.103) 


Dimuon 


174 


1.289 (0.941) 


1.018 (0.743) 


1.000 (0.730) 


0.996 (0.727) 


1.241 (0.906) 


Inclusive I 


194 


0.978 (0.711) 


0.961 (0.699) 


1.000 (0.727) 


1.029 (0.748) 


1.085 (0.789) 


Dimuon+I 


368 


1.126 (0.820) 


0.989 (0.720) 


1.000 (0.728) 


1.014 (0.738) 


1.160 (0.844) 


Inclusive II 


2097 


1.012 (1.120) 


1.014 (1.123) 


1.000 (1.107) 


1.011 (1.119) 


1.011 (1.119) 



Table II: Fit results using the LO dimuon calculation. We present the integrated strange-quark asymmetry [S ], and 
the cells display xVxoj an d the values in parentheses are the \ 2 1 DoF for each data subset. Note Xo is calibrated 

from the LO 5-fit. 



NLO Dimuon 


#pts 


B+ 


A 


B 


C 


B~ 


[S~] x 10 4 


- 


63.75 


23.92 


13.72 


12.81 


-18.59 


Entire Data 


2465 


1.033 (1.120) 


0.995 (1.079) 


0.994 (1.077) 


0.994 (1.077) 


1.025 (1.111) 


Dimuon 


174 


1.660 (1.212) 


0.960 (0.701) 


0.936 (0.683) 


0.937 (0.684) 


1.416 (1.034) 


Inclusive I 


194 


0.977 (0.710) 


0.974 (0.708) 


1.008 (0.733) 


1.017 (0.739) 


1.088 (0.791) 


Dimuon+I 


368 


1.301 (0.947) 


0.968 (0.705) 


0.974 (0.709) 


0.979 (0.713) 


1.245 (0.906) 


Inclusive II 


2097 


1.012 (1.120) 


1.009 (1.117) 


1.007 (1.115) 


1.005 (1.113) 


1.010 (1.118) 



Table III: Fit results using the NLO dimuon calculation. We present the integrated strange-quark asymmetry [S 
and the cells display x 2 /Xqi an d t ne values in parentheses are the \ 2 /DoF for each data subset. 
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